Precise calibration of cavity forward and reflected signals using low-level radio-frequency system* 
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Precise measurements of the cavity forward (Vf) and reflected signals (V,) are essential for characterizing 
other key parameters such as the cavity detuning and forward power. In practice, it is challenging to measure 
Vg and V; precisely because of crosstalk between the forward and reflected channels (e.g., coupling between the 
cavity reflected and forward signals in a directional coupler with limited directivity). For DESY, a method based 
on the cavity differential equation was proposed to precisely calibrate the actual V; and V,. In this study, we 
verified the validity and practicability of this approach for the Chinese ADS front-end demo superconducting 
linac (CAFe) facility at the Institute of Modern Physics and a compact energy recovery linac (CERL) test ma- 
chine at KEK. At the CAFe facility, we successfully calibrated the actual Vẹ signal using this method. The result 
demonstrated that the directivity of directional couplers might seriously affect the accuracy of Vs measurement. 
At the cERL facility, we calibrated the Lorentz force detuning (LFD) using the actual V;. Our study confirmed 
that the precise calibration of V; significantly improves the accuracy of the cavity LFD measurement. 
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I. INTRODUCTION 


The Chinese ADS front-end demo superconducting linac 
(CAFe) was constructed at the Institute of Modern Physics 
to develop the key technologies for a superconducting (SC) 
front-end linac. It was used to demonstrate the possibility of 
a 10 mA high-power continuous-wave (CW) proton beam for 
the China Initiative Accelerator-Driven System project [1-3]. 
The CAFe facility is a 162.5 MHz SC radio frequency (RF) 
facility operating in CW mode and consists of a normal con- 
ducting section and an SC section. As shown in Fig. 1, the 
normal conducting section is composed of an ion source, low- 
energy beam transport line (LEBT), RF quadrupole accelera- 
tor, and medium-energy beam transport line (MEBT). The SC 
section is composed of SC accelerating units, which consist 
of 23 SC half-wave resonator cavities and are assembled in 
four cryomodules (CM1 to CM4). 

The cavity forward and reflected signals (V¢ and V,) are 
typically used to control other key parameters such as the 
cavity voltage signal (V.), cavity detuning (A f), and forward 
power. For example, at the SC test facility (STF) at KEK [4], 
because there is no pickup antenna in the STF RF gun cav- 
ity, Vo is reconstructed by computing the superposition of Vp 
and V; (Ve = Ve + Ve) [5]. Note that in this paper, Ve, Vg, 
and V, refer to baseband signals, which can be represented by 
their amplitude and phase (e.g., Ve = |Ve|e7<""). The cavity 
pickup signal at the Paul Scherrer Institute was found to be 
corrupted by noise [6]. Therefore, Vs was recalculated using 
the actual V; and V,, and the noise was successfully removed. 
In addition, Vp measurement can affect the accuracy of the 
cavity dynamical detuning. For instance, at the compact en- 
ergy recovery linac (CERL) test facility at KEK, we compared 
the cavity Lorentz force detuning (LFD, A fLFp) calculated 
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using the cavity differential equation with that calibrated us- 
ing the square of the accelerating gradient. We found that the 
results of these two methods do not agree at the beginning 
of cavity filling. This discrepancy was attributed to the inac- 
curacy of Vp measurement, which was caused by the limited 
directivity of the RF directional coupler [7]. 


At CAFe, we measured the Ve, Vs, and V; signals and the 
low-level RF (LLRF) output in cavity CM3_5 (see Fig. 2). 
The location of the LLRF output signal is marked in Fig. 3. 
The cavity amplitude loop was operated in closed-loop mode, 
whereas the phase loop was operated in open-loop mode. As 
shown in Fig. 2, the cavity amplitude error was well compen- 
sated by the feedback loop. By contrast, because the phase 
loop was open, the cavity phase fluctuated greatly owing to 
perturbations such as microphonics and A fupp. Moreover, 
Fig. 2 shows that there is a significant difference between the 
V; phase and LLRF output phase. One possible reason for the 
phase difference is the nonlinearity of the solid-state amplifier 
(SSA). Another possible reason involves the inaccurate mea- 
surement of V; owing to crosstalk between the forward and 
reflected channels, for example, coupling between the cavity 
reflected and forward signals in the directional coupler. Later 
sections of this article demonstrate that the nonlinearity of the 
SSA is not the primary reason for this phase difference. 


Because of crosstalk between the measurement channels, 
the measured Vç signal is mixed with the V, signal, result- 
ing in inaccurate Vp measurement. To address this problem, 
Brandt proposed a method of eliminating the effect of V, on 
the V; signal [8]. We calibrated the actual value of Vp in 
CM3~_5 at CAFe and in an SC cavity at KEK-cERL using this 
method. The measurement results confirmed that this method 
is effective. 


The remainder of this paper is organized as follows. Sec- 
tion II briefly describes the CAFe LLRF system. Section II 
reviews the algorithms used to calibrate the actual values of 
V; and V,. In addition, the actual V; and V, values are esti- 
mated with higher precision. Section IV demonstrates the va- 
lidity of the method using the experimental results for CAFe 


and cERL. Section V presents a discussion of our future work. 
Finally, we summarize our study in Sec. VI. 


Il. LLRF SYSTEM 


At CAFe, the LLRF data acquisition system is used to mea- 
sure the amplitude, phase, and frequency and then use vari- 
ous signal processing methods to monitor and control the RF 
fields. Figure 3 shows a simplified illustration of the layout of 
the digital LLRF system at CAFe. The frequencies of the RF 
signals and the intermediate frequency (IF) signal are 162.5 
and 25 MHz, respectively. A field-programmable gate array 
(FPGA) module is used for real-time signal processing and 
high-speed data acquisition. RF signals such as Ve, V¢, and V; 
are first converted to IF signals. The IF signals are sampled 
at 100 MHz via a 16-bit analog-to-digital converter (ADC) 
and are fed to the FPGA. Then, the amplitude and phase of 
Ve are extracted from the IF signals. Subsequently, the ampli- 
tude and phase signals are compared with their setpoints, and 
their errors are calculated. These error signals are controlled 
by a proportional and integral (PI) controller. The controlled 
amplitude and phase signals are used to rebuild the IF signal. 
Next, the IF signal is up-converted, and the up-converted sig- 
nal is used to drive the SSA and ultimately power the cavity. 
An ARM chip is integrated into the FPGA. An experimental 
physics and industrial control system (EPICS) is installed on 
the ARM chip as an embedded Linux system for data acqui- 
sition. The amplitude and phase waveforms of Ve, Ve, and 
V, are saved to the hard drive of the CPU controller during 
measurement, and the measurement data are analyzed offline. 


Il. CALIBRATION OF REAL CAVITY FORWARD 
SIGNAL 


To calibrate the real V; and V, signals (see Fig. 2), we re- 
view the method of Brandt and present an example for cavity 
CM3_5 at CAFe. CM3_5 is mounted in the third cryomodule 
(CMs). It is the fifth cavity in CM3, and the loaded quality 
factor is 3.8 x 10°, as illustrated in Fig. 1. 

Figure 4 shows the V., V;, and V, signals measured by the 
LLRF system during the pulse conditioning of cavity CM3_5. 
According to these results, the cavity signal cannot be ob- 
tained by computing the measured signal between the V,* and 
V,* vectors if the two signals are uncalibrated. The actual cav- 
ity signal V; is a superposition of the actual V; and V; signals, 
and it can be reconstructed using the formula [8—10] 


Ve = Ve + Vi. (1) 
The actual V; and V, signals can be calibrated using the 


measured cavity forward and reflected signals according to 
[5, 8] 


aak o 


where V," and V* are the measured values. In this equation, 
a, b, c, and d are complex constants related to the actual and 
measured values. Then, Eq. (1) can be written as 


V = XVF + YV“, (3) 


where X = a + cand Y = b + d. The complex constants 
X and Y can be estimated by the multiple linear regression 
method. The specific calculation methods are as follows. 

The complex signals V, X, and Y can be written as V = 
R(V) + 91(V), X = Cy + jC2, and Y = C3 + 7C4, respec- 
tively. R(V), Cı, and C3 represent the real parts, whereas 
I(V), C2, and C4 represent the imaginary parts. Then, Eq. (3) 
can be written as 


R(Ve) + j1(Ve) = (C1 + JC2)[R(VE) + jI) 
+(C3 + jCa)[R(V,") + JIV). 
For a single pulse with sampled data points Vek, Viz, and 


Vz. where k runs from 1 to n, this method can be used to 


obtain a linear, overdetermined system of equations. Let V 
be the response variable; then 


y = [R(Vo1) IVen)". (5) 


In addition, let the measured values VF and V,* form the re- 
gressor matrix 


(4) 


R(Von) I(Ve1) yos 


RV) —I(Vä) RWI) —I(Vi) 
R(VŽ) —1(Ve) R(V3) —I(V3) 
IVE) ROVE) T(V) ROVA) 
2) ROZ) I(V3) RVŠ) 
T(Vin) ROR) IVa) ROn) 
Then the relation 
y=AC (7) 
is approximately satisfied for some parameter vector 
C=[01 Cz Cs Ga)”, (8) 


where C represents the values of interest. The best estimates 
of C1, C2, C3, and C4 can be determined using multiple linear 
regression. The complex constants X and Y can be obtained 
using these C values. The crosstalk elements b can be taken 
into account using the information from the RF decay time 
(after 3.5 ms in Fig. 4). During decay, the actual Vf should be 
0, and the factors a and b are then limited by another complex 
factor Z as 


Z=a/b=-—V,"/Ve (t> 3.5ms). (9) 


Then the complex factors can be expressed in terms of X, 
Y, and Z according to 


a=a 

b=a/Z (10) 
c=X-a 

d=Y —-a/Z. 
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Fig. 1. (Color online) Layout of CAFe facility. The half-wave SC cavities are mounted in four cryomodules (CM1 to CM4). Cavity CM3_5 


is marked by a red triangle. 
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Fig. 2. (Color online) Comparison of amplitude and phase of cavity voltage signal (V.), cavity forward signal (Vs), cavity reflected signal (V,), 
and LLRF output signal. The cavity amplitude loop was operated in closed-loop mode, whereas the phase loop was operated in open-loop 
mode. Figure shows that the phase of V; differs from that of the LLRF output. 


The final unknown complex factor a can be extracted nu- 
merically from the cavity differential equation assuming a 
constant cavity half-bandwidth wo.5, which can be estimated 
in terms of the field decay curve (t >3.5 ms). 

The cavity equation can be written in polar coordinates as 


d|Ve| 


di (11) 


+ |Velwo.5 = wo.5 : (2|Ve|)cos(0 — ¢). 


Consequently, the differential of |V.| can be computed using 
the formula 


y(a, tk) = wo.s [(2|Vs[)cos(@— 9) —|Vel], aD 


where 0 and ¢ represent the phases of V, and Vp, respectively. 


The variable y(a) represents the calibrated a Vel as a function 
d| Vel 


of the factor a according to Eq. (13). By contrast, ~q can 
be directly calibrated from the amplitude of V., which is inde- 
pendent of a. The difference between these two calculations 
should be minimized by optimizing a. Therefore, we denote 
the difference as à? (a). 
2 
ula) 


div 
MORDI ka - (13) 


tk 
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Fig. 3. Schematic of digital LLRF system at CAFe. The location of 
the LLRF output signal is marked. 
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Fig. 4. (Color online) Measured cavity voltage, forward, and re- 
flected signals in CM3_5 during RF commissioning. 


Note that a is a complex number with two components. A 
two-dimensional map can be obtained by sweeping the real 
and imaginary components of a, as shown in Fig. 5. Here, 
the color axis represents the value of \?(a). The X and Y 
axes represent the real and imaginary components of (a/ X), 
respectively. The optimal value of (a/X) is clearly 0.99 + 
0.072. To better illustrate this result, Figs. 6 and 7 compare 
the waveforms of V; and y(a), respectively, for a/X = 1 and 


a/X = 0.99 + 0.07%. The difference between alvel and y is 


minimized when the optimal value of a is used. 
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Fig. 5. (Color online) Scanning result of the parameter \?. The 
optimized factor a/X is 0.99 + 0.07i. The X and Y axes represent 
the real and imaginary parts of a/ X, respectively. 
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Fig. 6. (Color online) Comparison of cavity forward voltage V; for 
a/X = 1anda/X = 0.99 + 0.07%. The waveforms of XV and 
V; are also presented for comparison. 


IV. EXPERIMENTAL VERIFICATION 


To demonstrate the method described in the previous sec- 
tion, the V measurement data were processed and analyzed 
at the CAFe and cERL facilities. 
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Fig. 7. (Color online) Comparison of parameter y for a/X = 1 


and a/X = 0.99 + 0.07i. The quantity dl is also plotted as a 
reference. 


A. Experimental result at CAFe 


During the RF commissioning at CAFe, the LLRF output, 
SSA output, and V* and V,* signals of CM3_5 were down- 
converted and sent to the ADCs of the LLRF system. The 
positions of the four signals (points A, B, C, and D, respec- 
tively), are shown in Fig. 8. As mentioned in Sect. I, the 
phase of V,* differs from the LLRF output phase according to 
the measurement (see Figs. 2 and 9). In theory, the LLRF out- 
put phase and V; phase should be in agreement, whereas the 
phases at points A and C differ greatly, according to Fig. 9. 
As shown in Fig. 8, in the ideal case, the phases of the LLRF 
output (point A) and V; (point C) should be in agreement. The 
discrepancy has two possible causes: the nonlinearity of the 
SSA and crosstalk between the measurement channels (that 
is, the unwanted coupling of some components of V, with 
Vg). We analyze the effects of these two causes as follows. 

To eliminate the impact of the nonlinearity of the SSA, we 
measured the output characteristics of the SSA, as shown in 
Fig. 10. The amplitude and phase characteristics of the SSA 
output were measured using a network analyzer. The oper- 
ating point for Fig. 2 (and Fig. 9) is marked on Fig. 10 as 
P(x,y). The amplitude nonlinearity curve can be fitted using 
the polynomial formula [11] 


fala) = hns” + hnit +++ + he + ho. (14) 
The phase characteristic curves can be fitted by 
fo(x) = qm” + Om—1e™ 1 +- +r +q. (15) 


The amplitude (fa (|ViLRF|)) and phase (f(|ViLRF|)) of the 
SSA output can be calibrated using Eqs. (14) and (15), as 
shown in Fig. 9. 


Furthermore, fa (|VLLrRF]) and fo(|ViLRF|) can also be es- 
timated using the differential near the operating point. If the 
SSA is linear, its output amplitude increases linearly as the 
input amplitude increases, and the linear gain is K = y/a. In 
fact, the output amplitude of the SSA is often nonlinear with 
respect to the input amplitude; thus, the differential at point 
P is approximately Œ = Ay/Az, as illustrated in Fig. 9 
[11-13]. The amplitude variation of the SSA output signal 
depends on the value of G/K. For G/K > 1, the amplitude 
fluctuation of the LLRF output will be G/K times larger than 
that of the SSA. By contrast, for G/K < 1, the amplitude 
fluctuation of the LLRF output will be G/K times smaller 
than that of the SSA. We can see from Fig. 10 that the value 
of G/K is approximately 1.8; thus, a 1% fluctuation in the 
LLRF output amplitude near the operating point corresponds 
to a 1.8% fluctuation in the SSA output amplitude. It is easy 
to see that the ratio of the LLRF output amplitude and f4 is in 
good agreement with the value of G'/K in Fig. 9. In addition, 
a 1% fluctuation in the LLRF output amplitude near the oper- 
ating point will result in an error of 0.015° in the SSA output 
phase, according to the phase curves in Fig. 9. Noise from 
various sources also contributes to the error, which is difficult 
to determine from the LLRF output phase, as shown in Fig. 9. 
Thus far, we have estimated the amplitude and phase at point 
B (fa and fọ) successfully. Figure 9 shows that fa and fe dif- 
fer greatly in amplitude and phase (especially in phase) from 
Věž. Thus, we deduce that the nonlinearity of the SSA is not 
the main reason for the disagreement in phase by comparing 
the phase curves of fo(|ViLRF|) and fa. 


To determine whether the V, signal is coupled with the 
Ve signal, it is necessary to calculate the actual V; value of 
CM3_5. The values of the complex constants a and b must 
be known before this calculation can be performed. The ex- 
act value of a for CM3_5 was obtained by theoretical analysis 
and derivation in Sect. III. The exact value of b can then be 
calculated using Eq. (9). From the known a and b values, 
the actual cavity forward signal can be calibrated using the 
measured signals V,* and V,*. Figure 11 shows the actual Vę 
signal, measured Vf signal, and SSA output signal consider- 
ing the effects of the nonlinearity of the SSA. Panels (a), (b), 
and (c) represent the measurement results using the same cal- 
ibration factors (i.e., a and b) but for different dates. The red 
lines represent the amplitude and phase of the SSA output sig- 
nal considering the effect of the nonlinearity of the SSA. The 
blue lines show the calibrated actual forward signal, which 
represents Vp in the figure. The aqua lines represent the mea- 
sured forward signal, as in the previous section, which repre- 
sents Ve. Figure 11(a) shows good phase agreement between 
the blue and red lines. The LLRF output phase is in agree- 
ment with the phase of the V; signal in theory. This result 
confirms that the V, signal is coupled with the V; signal in the 
signal measurement. It is demonstrated that crosstalk in the 
measurement channels (e.g., crosstalk caused by the limited 
directivity of the directional coupler) is the major reason for 
the inaccurate measurement of V¢, and thus the phase discrep- 
ancy. 


To further confirm the results of V¢ signal calibration, we 
used the same calibration array (that is, the same values of 
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Fig. 8. Diagram of signal measurement. Points A and B represent the LLRF and SSA output signals, respectively. Points C and D represent 
the measured values of the cavity forward and reflected signals (V;* and V“), respectively. 
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Fig. 9. (Color online) Amplitude and phase of LLRF output signal, 
which corresponds to point A in Fig. 8 (green, VLLRF), SSA output 
signal considering the effect of SSA nonlinearity, which corresponds 
to point B in in Fig. 8 (red), and forward signal of measured corre- 
sponds to V (aqua). When the nonlinearity of the SSA is taken into 
account, the phases at point B (fo(|ViLRF|)) and point A (ZV;") still 
differ greatly. 


a, b, c, and d) to calibrate the measurement results for other 
cases. In these cases, the working point of the SSA is the 
same as the working point in Fig. 10. The blue and red lines 
show good agreement [see Fig. 11(b) and 11(c)] in all cases, 
indicating that the calibration array has high universality. 


B. Experimental result at KEK-cERL 


The cERL facility was constructed at KEK to accommo- 
date industrial applications of SC technology. It is a 1.3 GHz 
SC machine operated in CW mode [14-16]. At cERL, two 
nine-cell cavities (ML1 and ML2) with high Q, values of up 
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Fig. 10. (Color online) Schematic of amplitude and phase character- 
istics of SSA. The SSA and LLRF output amplitudes are normalized 
to the saturation point of the SSA. 


to 1x10" were installed in the main linac. 


At cERL, we observed that the LFD calculated using Vę 
differs from the theoretical value at the beginning of cavity 
filling in ML2. The discrepancy was assumed to be related to 
the inaccuracy of the Ve measurement due to signal crosstalk 
between the measurement channels. To verify the correctness 
of this assumption, the experimental data were re-analyzed. 
The cavity detuning (Aw) during the RF pulse without the 
beam can be estimated using the cavity differential equation 
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Fig. 11. (Color online) Amplitude and phase of SSA output signal considering the effects of the nonlinearity of the SSA (red), Viz (the actual 
forward signal, blue), Vi (the measured forward signal, aqua). Panels (a), (b) and (c) show the measurement results obtained using the same 
calibration factors a and b but for different dates. These cases show good phase consistency between the blue and red lines. 


[7, 8, 12] 


_ d _ wo.5(|2Vs|)sin(@ — ¢) 
dt [Ve] i 


Aw (16) 


where ¢ and 0 represent the phases of V. and Vç, respectively. 
The quantity wo.5 represents the cavity half-bandwidth. The 
detuning calculated using Eq. (16) includes other sources of 
detuning in addition to the LFD, although the LFD plays a 
major role in the overall detuning when the SC cavity is op- 
erated in pulse mode. 

Moreover, for a given pulse pattern, the LFD can be cali- 
brated using the LFD transfer function model K (s) [7, 17- 
19] as 


AWLFD = —2r K (s) E? 


acc? 


(17) 


where K(s) is the LFD transfer function, and E,,, is the 
accelerating gradient. The detuning at the beginning of the 
filling time can be estimated using the above equation [see 
A fmoaet in Fig. 12(b)]. 

According to Eq. (16), the amplitude and phase (|V;| and 
0, respectively) of Vç are required to calculate the cavity de- 
tuning Aw. When the crosstalk components b and c are not 
considered, the calibrated cavity forward signal Vẹ is equal to 
XV,* [Calil in Fig. 12(a)]. Under this assumption, the detun- 
ing A fcail is calibrated as shown in Fig. 12(b). Note that in 
Fig. 12(a), the amplitude of the Calil curve does not decrease 
to zero immediately when the RF power is switched off. As 


shown in Fig. 12(b), the waveforms of A fCail1 and A fmodel 
(red curve) are slightly different at the beginning of the filling 
time. 

Next, the crosstalk coefficients b and c in Eq. (2) are con- 
sidered, and the coefficients a, b, c, and d are recalculated us- 
ing the method described in Sec. II. Then, the actual cavity 
forward signal V; is calibrated by Vp = aV,* + bV,* [Cali2 in 
Fig. 12(a)]. As shown in Fig. 12(a), the amplitude of the Cali2 
curve decreases to zero immediately when the RF power is 
turned off. Using the waveform of Vç in the cali2 curve, the 
LFD (A feaiiz) is obtained, as shown in Fig. 12(b). As ex- 
pected, the trends of A fcaiig and A fmodel are in good agree- 
ment at the beginning of the filling time. 

From the analysis above, we conclude that the inaccurate 
measurement of V¢ can result in the miscalculation of the 
LED. In addition, the analysis demonstrates the validity of 
the method of Brandt. 


V. FUTURE WORK 


In the future, after estimating the values of a, b, c, and d 
in the offline state, we will implement Eq. (2) in the FPGA. 
The actual cavity forward signal Vẹ and the actual reflected 
signal V, are expected to be calibrated in real time. In the 
next stage, the calibrated V; can be used to calibrate other 
parameters such as the cavity detuning and beam current. 

In addition, the directivity of the directional coupler is re- 
lated to its input power (the output of the RF source); there- 
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Fig. 12. (Color online) (a) Two methods of calibrating the actual 
cavity forward signal. Green line, Calil: calibration of V; without 
considering signal crosstalk. Blue line, Cali2: calibration of V¢ con- 
sidering signal crosstalk. (b) Comparison of LFD calculated using 
Eq. (16) without considering signal crosstalk (green), using Eq. (16) 
considering signal crosstalk (blue), and using Eq. (17) (red). 


fore, the value of a/X may also depend on the location of the 
nonlinear operating point of the SSA. In the next step, we will 
attempt to determine the dependence of the parameter a/X 
on the operation point of the RF source. This research will 
be based on many measurements and calculations at CAFe at 
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VI. SUMMARY 
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